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THE DOPEX CODE - AN APPLICATION OF THE METHOD OF 
STEEPEST DESCENT TO LAMINATED -SHIELD -WEIGHT 
OPTIMIZATION WITH SEVERAL CONSTRAINTS 
by Gerald P. Lahti 
Lewis Research Center 

SUMMARY 

A two- or three- constraint, two-dimensional radiation shield weight optimization 
procedure and a computer program, DOPEX, is described. The DOPEX code uses the 
steepest descent method to alter a set of initial (input) thicknesses for a shield config- 
uration to achieve a minimum weight while simultaneously satisfying dose constraints. 
The code assumes an exponential dose- shield thickness relation with parameters spec- 
ified by the code user. The code also assumes that dose rates in each principal direc- 
tion are dependent only on shielding thicknesses in that direction. The code will pres- 
ently handle problems with one constraint (spherical geometry), or two or three 
constraints (right circular cylindrical geometry). Other geometries may be used if one 
subroutine is changed to suit the application. Code input instructions, a FORTRAN IV 
listing, and a sample problem are given. Typical computer time required to optimize 
a seven- layer shield is less than 0. 1 minute on an IBM 7094- n. 

INTRODUCTION 

One- dimensional shield weight optimization procedures have been established for 
radiation shield design (refs. 1 and 2). Optimization procedures such as these take 
some initial configuration specified by the user and alter the thickness of each layer to 
achieve a minimum weight while simultaneously satisfying a dose constraint. One 
method used to seek the minimum weight configuration from a starting configuration is 
that of steepest descent. 

Bernick (ref. 3) developed an algorithm to optimize one- dimensional shields utiliz- 
ing the steepest descent method which was programmed as the OPEX code (ref. 4). 



OPEX was later reprogrammed as the OPEX n code (ref. 5). This report presents the 
DOPEX code, an extension of this type of procedure to two dimensions with two or three 
constraints. The statement of the optimization algorithm and the derivation of the ex- 
pression for the steepest descent vector were made in reference 6 for the general multi- 
dimensional, multiconstraint case. The present procedure considers only the special 
case where doses in each principal direction are dependent only upon the shielding thick- 
ness in that direction. This particular dose-thickness model is valid for many config- 
urations of interest. 

The geometry and thicknesses of material determine the weight and the derivative 
of weight with respect to thickness. With first derivatives of both weight and dose with 
respect to thickness, the optimization procedure alters the base configuration to obtain 
a set of thicknesses corresponding to a minimum weight configuration while meeting the 
dose constraints. This code does not select materials but merely alters initial config- 
urations. DOPEX may eliminate layers but it never adds layers. The accuracy of the 
optimization procedure depends on the ability of the dose- thickness relation to accurate- 
ly represent the actual dose and changes of dose with thickness in the real geometry. 
Although the present model has been reasonably successful, a final detailed proof cal- 
culation is necessary to confirm dose predictions of the optimization code for the re- 
sultant optimized geometry. It is realized that some shield shaping (e.g. , rounding of 
the corners) may result in further weight reduction over that predicted by an optimiza- 
tion. But such a procedure is beyond the scope of this report. 

The two-dimensional shield- weight optimization code DOPEX* described here is 
designed for one, two, or three constraints. The subroutine WEIGHT included in this 
report is presently programmed for one- constraint spherical geometry (same as in 
ref. 5) and two- and three- constraint right circular cylinders. Any different geometry 
that is being considered requires writing a new subroutine WEIGHT for that geometry. 

In this report, a description of how to obtain the necessary input data for DOPEX is 
given. Data input instructions, a FORTRAN IV code listing, and a sample problem op- 
timizing a three- constraint seven- layer shield of tungsten and lithium hydride for a 
space power reactor are given. 


DOSE-THICKNESS RELATION 

The total radiation dose rate D m at some reference point in space is defined, for 
purposes of the optimization procedure, to be 

*In the evolution of this code from the original OPEX code, one version was called 
D'OPEX standing for double- constraint OPEX. Although the code was extended to three 
constraints, the name DOPEX was retained. 
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IMAX m 



(1) 


where 

tVi 4-u 

D im i 1 component of m in total dose rate (e. g. , dose due to capture gammas 
from first shield layer or dose due to inelastic gammas from last shield 
layer, etc. ) 

D m total dose rate in m^ h direction (i. e. , at m th constraint point) 

J-U 

IMAX^ number of components of dose in m L direction F 

Each dose component D im is further assumed to be of the form 


NREG, 


D. = C- exp 
lm 1m ^ 


I 




1=1 


u.. t- 
j M ijm jm | 


( 2 ) 


where 

C im 

NRECm 

t- 
]m 

M ijm 

Subdividing the total dose rate into its components and evaluating coefficients M^j m 
separately for each component dose rate results in a total dose model having a set of 
coefficients which is relatively insensitive to changes in thicknesses of layers. (See, 
e. g. , the approach used in ref. 1. ) Other forms of the dose- thickness relation could be 
utilized but this would require reprogramming the subroutine which calculates dose, and 
dose gradients. The parameters in this empirical expression are assumed to be ob- 
tained by fitting to accurate, detailed radiation transport calculations of dose rate for a 
given base configuration and perturbations of that configuration. 

Notice from equation (2) that the nr dose is alterable only by thicknesses in that 
direction. This assumption was made in the final derivation of the equations in refer- 
ence 6 because it may be virtually impossible to write a dose equation such as equa- 
tion (2) involving n-dimensional shields. The user is cautioned to be aware of the clas- 
sical shielding booby traps of very asymmetric shields which may be generated by this 
code (or any other such code) which uses a 'Tine- of- sight" type of representation of the 


fitted parameter 

number of regions present in m direction 

il x i_ 

thickness of j region in m l direction 

"attenuation coefficient" which describes effect of change in thickness tj m 
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dose- thickness model. Although the doses in the various directions are decoupled from 
one another, the derivatives of weight, dose, and products of these are not. Note that 
there may be more than one contribution to the total dose from a given layer. Other 
layers may contribute negligible secondary gamma dose and their dose contribution is 
omitted. Generally because of differences in formation rate and gamma attenuation, 
a given high-Z shield region will have one dose component due to capture gammas and 
a second due to inelastic gammas. The core will have primary neutrons and primary 
gammas. A hydrogeneous layer may have negligible dose. 

For example, consider the cylindrical reactor- shield described in sketch (a) and 
the following table with a dose constraint at points P^ and Pg: 


• p 


1 



Region j 
(both di- 
rections) 

Description 

Radiation 
dose com- 
ponent i 

Description 

1 

Reactor core 

1 

Core neutron 



2 

Core gamma 

2 

Reflector 

3 

Capture gamma 



4 

Inelastic gamma 

3 

Lithium hydride 

— 

Negligible dose 

4 

Tungsten shield 

5 

Capture gamma 



6 

Inelastic gamma 

5 

Lithium hydride 

— 

Negligible dose 
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In this example, jUg 4 > M 3 45 and M 4 4 lor directions 1 and 2 represent the atten- 
uation of gamma rays by the tungsten layer for the various gamma sources in the core 
and layers between the core and the tungsten. The coefficient ^ represents the 
attenuation by the tungsten of neutron dose due to neutrons born in the core. In contrast, 
ju^ 2 represents the attenuation by lithium hydride in region 3 of neutrons which give 
rise to capture gamma sources in the tungsten (region 4). 

Coefficients, such as ju^ describe the effect of the thickness of a source region 
(in this case tungsten) on the dose from its own source. Equation (2) becomes inade- 
quate particularly when the thickness of this layer goes to zero. DOPEX does, however, 
set C im equal to zero if a particular layer is eliminated by the code. At this point, 
however, the attenuation coefficients u, • should be recalculated for the new config- 

uration. Regions of origin of the various dose components are required data input to 
the code to facilitate this operation. 


FITTING PARAMETERS TO DOSE-THICKNESS RELATION 

The coefficients of equation (2) are obtained as precisely as possible by performing 
a series of transport calculations. A starting, base configuration is selected, prefer- 
ably as close to an optimum configuration as possible. With the transport calculation, 
such as a discrete ordinates calculation, the individual dose components D im are eval- 
uated. Each layer of the base configuration to be altered by the optimization code is 
then systematically increased by a nominal amount (say, 1 cm) and the dose components 
are reevaluated with the transport code. From equation (2), then, for each m and i 


^m^lm’ ^m’ * ' •» ^km’ • • •» t NREG m^ _ C im exp 


D im^lm’ ^2m’ ' ‘ ‘ ’ Han + at km> • • •> HiREG " C im exp 

ex P<-"ikm A W 

where NREG is the number of layers in the problem. Solving the previous pair of 
equations for results in 
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^ikm 


At km 


In 


D im ^lm’ t 2m’ * 


> ^ 


km’ 


D im ^lm’ ^m’ 


^km + A ^km’ 


t NREG 
• *» t NREG 


The coefficients jUjj m are determined for all i dose components, all j thicknesses, 
and all m directions in this manner. 

The basic configuration data, that is, the set of thicknesses t. , dose components 
D im , and "attenuation coefficients” constitute the required input data. For the 

base configuration, then, the coefficient C im is calculated by the code DOPEX from 
equation (2); that is, 



WEIGHT-MINIMIZATION PROCEDURE 

The procedure for obtaining the minimum weight configuration by the method of 
steepest descent is presented in this section. The equations are from reference 6. 

The mathematical problem to be solved is that of minimizing the weight w, a func- 
tion of all thicknesses tj m , while constraining the total doses D m to some particular 
values; that is, 

Minimize w(t u , t 21 , t 31 , . . ., t NREG 1? t J2 , t 22 , . . ., t^gg 2 , . . ., t^Q M ) 

with the following constraints: 

(1) D m = D m for m = 1, 2, . . . , M 

<*> *jm 2 0 

(-3) t^ = constant for any desired values of l 

where M is the total number of constraints and D° is the value of the dose constraint 
in the m direction. 

Constraint (2), non-negative shielding thickness, ensures a physically meaningful 
solution. The optional constraint (3), fixed thickness, is useful if it is desired that 
some thicknesses be kept from changing during the course of the calculation (e.g. , the 
reactor core and reflector thicknesses). Constraint (3) is necessary for the spherical 
and cylindrical geometry programmed into DOPEX to prevent the trivial case of reduc- 
ing the reactor core size to zero. 
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An n-dimensional Euclidean vector space with Cartesian coordinates t^, t 2 p 
. . . , p 1 22 , • • • , %reg M is defined - The following vectors are de- 

fined on this space : 

t£ ^ll’ *21’ • ' •» t NREGl» fc 12’ *22’ ' ' •» ^£0 2’- • * ^ 


— _ / 3w 3w 
g = 


3w 


3w 3w ' 


l3tll 3t 21 


> 7 

3t NREG 1 3t 12 3t 22 


at 


3w 

NREG My 


a =/ 9D 3D 

m 


3D 


3t lm 3t 2m 


at. 


NREG my 


The notation t = (t^, t 2 j, . . tjqREG means tjXj + t«jX 2 + ’ • • w h ere X: 

are unit vectors in the i directions. Vectors g and a represent the gradient of 
weight and dose, respectively. The components of g are evaluated from analytic ex- 
pressions of weight as a function of thickness and depend on geometry. The components 
of a are evaluated from the partial derivatives of (1) and (2), namely, 


3D, 


m 


3t km 


IMAX IMAX 

m m 

\ ' 3D. \ ' 3 

\ im _ \ 3 

Z . — J St km / J 3t km 


C. 

im 


i=l 


i=l 


exp 


NREG 

£ 


^ijm^mJ 


IMAX„ 


H 

i=l 




NREG 


^ikm> C im ex Pl- 


j=l 


u.. t- . 
ijm jm/ 


IMAX.. 




^ikm^im 


i=l 


The unit vector u (see ref. 6 for derivation) 

h • s 


• g 


-g + Z “ + — — a 2 + 


u = ■ 


a l • a l 


a 2 ’ a 2 


_ _ ( a l • g) 2 (a 2 • g) 2 

g • g — 


a l ' a l 


a 2 ‘ a 2 


+ a 
a M * a M 


M 


~\2 


( a M • g ) 


a M * a M 


1/2 
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points in the direction of greatest weight decrease (steepest descent) along a hyperplane 
tangent to the hypersurface described by the equation 

D m (t ) = constant = for all m = 1, 2, . . . , M 

Components of u represent increments of thickness to be added to each tj m to 
approach the minimum weight criterion. 

The optimization code proceeds as follows: 

(1) A fraction f of u is added to the thickness t. The fraction f is an input 
parameter. (A value of f = 1. 0 has given satisfactory results since |u| = 1. ) 

(2) The new set of thicknesses generally does not return the correct dose constraint 
so a first-order correction is applied to each t to return to the dose constraint. That 
is requiring 



= D m (t + dt) = D (t ) + VD (t ) • dt 
m m ' m ' 


one obtains (ref. 6) 


dt *[ D m- D mffl]:r 


m 


m 


m 


Finally 


Steps (1) and (2) are repeated until the relative change in weight from one iteration 
to the next is less than some prescribed value. 

The code output includes final thicknesses and individual dose components as calcu- 
lated from the dose-thickness relation for each iteration. It is incumbent on the user to 
make a final detailed proof calculation to verify the results of the prediction of the opti- 
mization code. Experiences have indicated that if input coefficients are determined 
accurately, final predictions of the optimization code are quite good provided the config- 
uration is not radically changed. If the configuration is changed severely, a recalcula- 
tion of coefficients is in order. 


THE DOPEX CODE 

In this section, the details of the DOPEX code are presented. Included here are: 
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(1) General comments 

(2) A flow chart for data input 

(3) A sample problem and sample problem output 

A complete FORTRAN IV listing appears in the appendix A. 

The subroutines and their function in the DOPEX code are as follows: 

(1) Main program - Calculates the steepest descent vector and alters shield thick- 

nesses; calls subroutines INPUT, INIT, DOZE, WEIGHT, and CLEAR 

(2) Subroutine INPUT - Reads input data 

(3) Subroutine INIT - If input thicknesses do not result in the same dose rate as 

that required by the constraint, INIT is called to calculate a set of thicknesses 
which does meet the constraint so that the optimization procedure can begin. 
This calculation is carried out without regard to weight minimization. Calls 
subroutines DOZE, WEIGHT, and CLEAR. 

(4) Subroutine DOZE - Calculates dose rates and derivatives of dose with respect 

to thickness and products of dose and weight derivatives. 

(5) Subroutine WEIGHT - Calculates the total weight and derivatives of weight with 

respect to thickness. 

(6) Subroutine CLEAR - If a shield layer is eliminated in the course of the calcu- 

lation, subroutine CLEAR is called to set that layer thickness to zero for the 
remainder of the calculation, set the constraint flag so that it is maintained at 
zero thickness, and zero out a dose component if it originated in the removed 
layer. 

The DOPEX code as presented here includes the following geometry- dose- constraint 
situations : 

(1) Spherical geometry, one dose constraint (this is the same as that presented in 

ref. 5) 

(2) Cylindrical geometry, two dose constraints 

(3) Cylindrical geometry, three dose constraints 

These geometries are illustrated in figure 1. Different geometries will require a re- 
coding of subroutine weight to suit the application. 

The cylindrical geometry (see fig. 1(b) and (c)) assumes the first dose constraint to 
be in the radial direction and the second and third to be on the ±z-axis as shown. Be- 
cause of the way in which subroutine WEIGHT is coded, the number of layers in each 
direction must be the same as in the other directions. Furthermore, the density of the 

i.L i.1. 

i layer in the first direction must be the same as that of the i layer in the second 
direction, and so on. If one has more shield layers in one direction than another, then 
one must artificially add layers such that there are the same number of contiguous re- 
gions in all directions. These added layers are fixed at zero thickness. Such a con- 
figuration is shown in figure 2. The present geometry will work with two or three 
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(a) One-constraint, spherical. 



(b) Two-constraint, cylindrical (four regions illustrated in each tc) Three-constraint, cylindrical (four regions i I lustrated), 

direction). 


Figure 1. - DOPEX geometry. 


r, constraint r Regions 2 and 3 have zero thickness in 
direction 1 / radial direction 
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Figure 2. - Cylindrical geometry with two constraints, illustrating region des- 
ignation to handle case of unequal nun, her of shield layers in different di- 
rections. 



constraints. The required geometry data for input consists of a set of thicknesses for 
each region and direction. Radii are calculated internally. A thickness must be spec- 
ified for each region. 

The code is constructed so that if one requires a new geometry, only a new sub- 
routine WEIGHT is required. This subroutine calculates the system weight and all 
weight derivatives. 

The code can presently handle up to 25 regions and dose components. Running 
times are typically less than 0. 1 minute on the IBM 7094-11 for a 12- region, 16 dose 
component problem . 


Flow Chart for Data Input 


The following is a flow chart for data input where the symbol 
read a card or set of cards as described: 



means 



Title card 

FORMAT (12A6) 


f 

Any alphanumeric information in card columns 1 to 72 

L 2 _ 

Control card 

FORMAT (215, 4E10. 4) 


Card 

column 

Variable 



1 to 5 

ICON 

Number of directions (constraints) in this problem: 
ICON = 1, 2, or 3 


6 to 10 

MAX 

Maximum number of iterations allowed 


11 to 20 

EPS 

Convergence criterion for weight 


21 to 31 

EPSD 

Convergence criterion for initial dose (EPS, EPSD 
typically 0. 001) 


31 to 40 

C0N 

Fractional step size on u (C0N ~ 1.0) 

1 

41 to 50 

CA 

Fractional step size for initial dose (CA ~ 1. 0) 

/ 

3 

Control card 

FORMAT (215, E10. 4) 


Card Variable 
column 

1 to 5 NREG Number of regions in first direction (<25) 
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6 to 10 IMAX 


Number of dose items in dose equation in first direc- 
tion (^25) 


11 to 20 DSTAR Value of dose constraint in first direction (used only 

if DSTAR >0.0) 



[ 3 tc 

l 


T(J) I FORMAT (7E1 0.4) 

it. 

Thickness of j Ln region in first direction; NREG values required. 

RH0(J) FORMAT (7E1 0.4) 

Density of j** 1 region in first direction; NREG values required. 

NB(J) thickness constraint flag FORMAT (2511) 

iv 

NB(J) = 0 constrains j region in first direction to constant thick- 
ness 

iL 

NB(J) = 1 allow j Ln region in first direction to change; NREG values 
required. 

Mjj FORMAT (7E1 0.4) 

For the first direction and for each region j, read new card (or set of 
cards) with IMAX values of ju^; NREG cards (or sets of cards) 
required. 


D(i) 

XL 

For the first direction, D(i) is i- dose component corresponding to 
initial geometry; IMAX values required. 

NIJ(I) FORMAT (1415) 

NIJ(I) is region number in first direction physically associated with 
i th dose term. This cross referencing is necessary to zero the 
correct C^ m in the event the region is diminished to zero thick- 
ness by optimization process. A region maybe identified more than 
once. IMAX values required. 


Read cards 3 to 9 containing information pertaining to second direc- 
tion. (Required if ICON >2. ) 


12 











Read cards 3 to 9 containing information pertaining to third direction. 
(Required if ICON = 3. ) 


Sample Problem 

A sample problem consisting of a reactor, molybdenum reflector, and a shield con- 
sisting of seven alternating layers of lithium hydride and tungsten is illustrated in fig- 
ure 3. To illustrate the use of this code, the problem is run for the three constraint, 


To dose constraint 1 
(2 mrem/hr) 


To dose constraint 3 
(20 mrem/hr) 



To dose constraint 2 
(2 mrem/hr) 


cylindrical geometry. Region descriptions, densities, and initial thicknesses (guessed) 
are listed in table I. The configuration is to satisfy a dose constraint of 2 millirem per 
hour at a point 20 meters distant in the radial (first) direction, 2 millirem per hour 
20 meters distant in the +z (second) direction, and 20 millirem per hour 20 meters dis- 
tant in the -z (third) direction. A series of spherical discrete ordinates calculations 
for both neutrons and gamma rays were made to calculate doses from each source re- 
gion for the initial set of thicknesses; perturbations were made to determine the atten- 
uation coefficients The results of these calculations are listed in tables II 
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TABLE I. - REGION DESCRIPTION 


Region j 

Description 

Density, 

g/cm 3 

Thickness, 

cm 

1 

Reactor core 

9.957 

26. 0 (radius), 




52. 0 (height) 

2 

Plenum 

8.647 

2. 50 

3 

Pressure vessel 

16.763 

.60 

4 

Molybdenum reflector 

9. 234 

11.00 

5 

Lithium hydride 

. 75 

a 17. 90 

6 

Tungsten 

19. 3 

a 7. 00 

7 

Lithium hydride 

. 75 

a 14. 00 

8 

Tungsten 

19.3 

a 5. 00 

9 

Lithium hydride 

.75 

a 10. 00 

10 

Tungsten 

19.3 

a 3. 50 


Lithium hydride 

. 75 

a 59. 50 


a Initial guess. 


TABLE II. - DOSE COMPONENTS CALCULATED 
WITH INITIAL 3 THICKNESSES 


Dose component 

Value of dose component, 
D i> 

mrem/hr 

1 Neutron 

0.02430 

2 Core gamma 

. 00303 

3 Plenum, pressure vessel capture gamma 

.00196 

4 Plenum, pressure vessel inelastic gamma 

. 00220 

5 Reflector capture gamma 

.204 

6 Reflector inelastic gamma 

. 00504 

7 Region 6 tungsten capture gamma 

.0921 

8 Region 6 tungsten inelastic gamma 

. 00974 

9 Region 8 tungsten capture gamma 

.0988 

10 Region 8 tungsten inelastic gamma 

.0278 

11 Region 10 tungsten capture gamma 

.201 

12 Region 10 tungsten inelastic gamma 

. 0947 

Total 

0.7647 


a From spherical calculations made with guessed shield thicknesses 
listed in table I. Assumed to apply in all directions of this cylin- 
drical optimization calculation. 
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and HI. Because the core radius, plenum, pressure vessel, and reflector thicknesses 
are to be constrained in this calculation, a value of /i.j m equal 0. 0 is assigned to these 
regions. The data calculated for the cylindrical case was determined from spherical 
discrete ordinates calculations and assumed to hold true for this cylindrical case; this 
may not, of course, be generally true. 

The computer output for this sample problem is listed in appendix B. The output 
consists, first, of a listing of all input information, followed by the value of the dose 
constraint (DES), the value of the calculated dose for the initial configuration (DOS), and 


TABLE IV. - FINAL SHIELD THICKNESS PREDICTIONS 
FOR OPTIMIZED CYLINDRICAL SHIELD 


Region 

Description 

Shield thickness, 

cm 

Radial 
direction 
(2 mrem/hr) 

+z 

direction 
(2 mrem/hr) 

- z 

direction 
(20 mrem/hr) 

5 

Lithium hydride 

20.37 

17. 19 

14.72 

6 

Tungsten 

9.48 

7.65 

6. 11 

7 

Lithium hydride 

9.76 

13. 12 

9.97 

8 

Tungsten 

2.20 

4. 14 

3.72 

9 

Lithium hydride 

11.28 

9.32 

6.26 

10 

Tungsten 

3.24 

2.80 

2.09 

11 

Lithium hydride 

41.35 

54. 39 

53.73 


the weight (WT) in grams of the initial configuration. If DES and DOS do not agree to 
within the parameter EPSD, a new set of thicknesses is calculated and printed out; this 
new set of thicknesses satisfies the dose constraint. The list of DES, DOS, and WT is 
followed by the values of the dose components D(i) and the thicknesses of each region 
T(i). For multidimensional calculations, the aforementioned blocks of data are repeated 
for each of the respective dose constraint directions as shown in figure 1. Following 
this, the results of each DOPEX iteration are listed. The data listed for each iteration 
are: 

IT, the iteration count number 

WT, the system weight (in units consistent with input data) 

DOS, the calculated dose rates in each of the constraint directions 
This is followed by the individual dose components D(i) and the region thicknesses T(i) 
for each of the respective constraint directions. The sample problem output shown in 
appendix B has been abbreviated to show the results of the first five and the last five 
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TABLE V. - DOSE COMPONENTS CALCULATED WITH OPTIMIZED THICKNESSES 


Dose component 

Value of dose component, 

D^, mrem/hr 

Radial 

direction 

+z direction 

-z direction 

1 Neutron 

0.2521 

0. 0735 

0. 4705 

2 Core gamma 

.0074 

. 0073 

. 0702 

3 Plenum, pressure vessel capture gamma 

.0046 

.0046 

. 0430 

4 Plenum, pressure vessel inelastic gamma 

. 0052 

. 0052 

. 0482 

5 Reflector capture gamma 

. 4418 

.4860 

4.9452 

6 Reflector inelastic gamma 

.0120 

. 0122 

. 1234 

7 Region 6 tungsten capture gamma 

.3562 

.3994 

4. 1715 

8 Region 6 tungsten inelastic gamma 

. 0798 

.0442 

.3090 

9 Region 8 tungsten capture gamma 

.2601 

.2801 

2. 8275 

10 Region 8 tungsten inelastic gamma 

.0629 

. 0767 

.7146 

11 Region 10 tungsten capture gamma 

.3298 

.4253 

4. 6897 

12 Region 10 tungsten inelastic gamma 

. 1880 

. 1856 

1. 5871 

Total 

2.0 

2.0 

20.0 


DOPEX iterations to avoid repetition. The final thicknesses and values of the dose com- 
ponents are listed in tables IV and V, respectively. The cylindrical problem here 
7 

weighed 5. 4X10 grams (119 000 lb) at the beginning of the optimization; the final weight 
was 4. 33X10 7 grams (95 000 lb). 


CONCLUDING REMARKS 

A two-dimensional, two- or three- constraint radiation shield weight optimization 
code DOPEX has been developed. The code uses the steepest descent method to seek 
out a minimum. An empirical dose- thickness model is used by DOPEX; the user sup- 
plies appropriately fitted coefficients. The model assumes that doses in each principal 
direction are dependent upon thicknesses only in that direction. Such a model is valid 
for many cases of interest. Typical problems require about 0. 1 minute of running time 
on an IBM 7094- n computer. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 6, 1972, 

112-27. 
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APPENDIX A 


FORTRAN IV LISTING 

The IBM 7094-n FORTRAN IV listing for the DOPEX code is presented in this ap- 
pendix. Subroutine WEIGHT is coded for the spherical and cylindrical geometries 
shown in figure 1. 


$IBFTC DOPEX C 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 


***D-OPEX *«*** A ONE-, TWO-, OR THREE CONSTRAINT RADIATION 
SHIELDING OPTIMIZATION CODE 

THIS CODE IS BASED ON THE ONE-C ONSTRA I NT CODE, 
OPEX-II, A STEEPEST-DESCENT-METHOD PROCEDURE. 
THIS PROCEDURE WAS ORIGINALLY DEVELOPED FOR 
ONE CONSTRAINT BY R. L . BERN I CK OF ATOMICS INTL. 

G.P. LAHTI NASA-LEWIS RESEARCH CENTER 

COMMON ICON, MAX, EPS, EPSD, CON, CA, WT, GG, IT, IC, 

1 DSTARI3), NREGI3), T(25,3J, RHO(25,3), NB(25,3), 

2 IMAXI3), 0(25,3). E MU ( 2 5, 2 5 , 3 ) , NIJ(25,3>, DUSE ( 3 ) , C(25,3>, 

3 A ( 25, 3 ) , A A ( 3 ) , AGO), U(25,3>, G(25,3) 


DIMENSION ALPHA ( 12 ) 

7777 READ ( 5 , 500 ) ALPHA 
WRI TE( 6,600) ALPHA 

RE AD (5, 505 ) I CON , MAX , EPS, EP SD, CON , C A 
WRITE( 6,605) I CON, MAX, EPS, EPSD, CON, C A 
DO 5 I C=l, ICON 

CALL I NPUT(NREG( IC) , I MAX ( IC ) ,OSTAR( IC ) , T(1,IC), RHOtl.IC), 
i NB(1,IC), EMU (1,1, I C ) , 0(1,10, NIJ(1,IC), C(1,IC)) 

5 CONTINUE 
C 

CALL WEIGHT 
WRI TE ( 6, 607 ) WT 


FOR EACH OF THE IC DIRECTIONS, CALCULATE DOSE AND NORMALIZE TO 
THE CONSTRAINT , IF NECESSARY 
DO 10 I C= 1 » I CON 

CALL DOZE(NREG( IC), IMAX(IC), T { 1 . I C ) . EMU(1,1,IC), NB(1,IC), 

1 CC.IC), 0(1,10, A ( 1 , 1 C ) , DOSt(IC), AA(IC), AG(IC), 6(1,10) 
IFIOSTARI IC) .LE. C.O) GO TO 9 

CALL INIT(NREG( IC), IMAX(IC), T(l.IC), EMU(1,1,IC), NB(1,IC), 

1 C(.,IC), D ( 1 , IC ) , A(l.IC), DOSE ( I C ) , AA(IC), AG(IC), G(l.IC), 

2 DSTAR(IC) ) 

GO TO 10 

9 DSTAR( IC)=DOSE( IC) 

10 CONTINUE 
C UPDATE ALL DERIVATIVES 

CALL WEIGHT 
DO : 5 I C= 1 , ICON 

CALL QOZE(NREG( IC), IMAX(IC), T(1,IC), EMUd.l.IC), NB(1,IC), 

1 C ( ^ , I C ) , 0(1,10, A ( 1 , I C ) , DOSE ( I C I , A A ( I C ) , AG ( I C ) , Gli.IC)) 
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15 CONTINUE 

C **«* THE ITERATION LOOP 
00 i 03 IT =i, MAX 
W = WT 

DO 25 IC = 1 , ICON 
SQ=GG 

DO 2D I K= 1 » I CON 
20 SQ=SQ-AG( IK)**2/AA( IX ) 

SQ=SQRT(SQ) 

NR=NREG( IC ) 

00 25 I =1, NR 

I F ( MB I I . I C ) .EG. 0) GO TO 25 
C INCREMENT EACH THICKNESS 

V= AG ( IC ) * A { It I C ) / A A ( IC ) -GII.IC) 

V=V/ S3 
U( I , IC )=V 

T( I , i: I =T ( I, I C I + CON*V 

IF(Td.IC) .LE. 0.0) CALL CL EAR ( I , I MAX ( I C ) ) 

25 CONTINUE 

DO lO I C= 1 « I CON 
NR=NRE G ( I C ) 

CALL OOZE (NREG( IC) , IMAX(IC), T(1,IC), EMU(l,i,IC), Nbll.IC), 

1 ClltlOt 0(1, IC), A ( 1 , I C ) , OOSE(IC), AA(IC), AG(IC), GI1.1C)) 

C RESTORE DOSE CONSTRAINT 

DO 30 1=1, NR 

I F ( NB( I , I C ) .EQ. 0) GO TO 30 
CONST = ( DS T AR ( I C ) -DOS £ ( IC ))/AA( IC) 

Til ,IC )=T( I , IC)+ CON S T*A ( I , I C ) 

1 F ( T ( I , IC ) .LE. 0.0) CALL CL E AR ( I , I MAX ( I C ) ) 

30 CONTINUE 

C UPDATE ALL DERIVATIVES 

CALL HEIGHT 
DO 35 I C= 1 * I CUN 

CALL DOZE ( NREGI IC ) , IMAX(IC), TCI, I C ) , EMU(1,1,IC). NBII.IC), 

1 C ( L , I C ) , D ( 1 , 1 C ) , A ( 1 , 1 C ) , DOSE(IC), AA(IC), AG(IC), G(l.IC)) 

35 CONTINUE 

C WRITE OUT RESULTS OF THIS ITERATION 

WRITE! 6,610) IT, WT, (OOSEIIC), IC=1,IC0N) 

DO A 5 I C= 1 « I CON 
I M= I MA X ( I C ) 

NR=NRE G ( IC ) 

WRITE(6,615) (0(1,10, 1=1, IM) 

45 WRITE! 6,620) (T(I,IC)» 1=1, NR) 

I F { A3SIWT-W)/ WT -EPS) 110, 110, 100 
100 CONTINUE 
110 GO TO 7777 
500 FORMAT (12A6) 

505 FORMAT (215, 4E10.4) 

600 FORMAT (1H1, 12A6) 

605 FORMAT { 9H0 ICON = 13/ 9H MAX = 13/ 

1 / 9 H EPS = 1PC12.4/ 9H EPSD = 1PE12.4/ 9H CON = 

2 1PEL 2,4/ 9H CA = 1PE12.4) 

607 FORMAT ( 46HU THE WEIGHT OF THE STARTING CONFIGURATION IS 1PE12.4) 
610 FORMAT ( 1H0/7H0 IT = I3,7X,5HWT = 1PE12.4.7X, 

X 7X, SHOOS = 1P3E12.4) 

615 FORMAT ( 7H0D ( I ) - 1 P 9E 1 2. 4/ { 7X, 1P9E12.4) ) 

620 FORMAT ( 7H0T ( I ) = 1P9E 1 2 . 4/ ( 7X , 1P9E 12. 4 ) ) 

END 
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$ I BF TC IN PUTT 


SUBROUTINE INPUT ( NR EG , I MAX , DSTAR , T , RHO , NB , EMU , D , N I J , C ) 

01 Mr Ns ION T(i), RHO(l), NB ( 1 ) , EMU(25,25), 0(1 ), NIJ(l), C(l) 
common ICON, MAX, EPS, EPSD, CON, CA, WT, GG, IT, IC 

READ ( 5 , 400 ) NREG, IMAX, DSTAR 
WRITE! 6,500) IC, NREG, IMAX, DSTAR 
READ (5,401) ( T ( I ) , 1=1, NREG) 

READ ( 5 » 4C 1 ) (RHO(I), 1=1, NREG) 

READ (5,403) ( NS ( I ) , 1=1, NREG) 

WRITE (6,510) ( I, T ( I ) , RHO(I),NB(I), 1 = 1, NREG) 

WRITE! 5,550) 

DO 7 J = 1 , NREG 

READ (5,401) (EMU(I,J), 1=1, IMAX) 

7 WRITE (6,570) J,(EMU(I,J), 1=1, IMAX) 

READ (5,401) ( D ( I ) , 1=1, IMAX) 

READ ( 5 , 402 ) (NIJ(I), 1=1, IMAX) 

C NIJ(I) IS THE REGION NUMBER ASSOCIATED WITH THE ITH DOSE TERM 

D0SE=3 • 0 
DO '■> 1=1, IMAX 
DOS_ = OOSE+D( I ) 

B8 = 3.0 
DO 3 J = 1 , NREG 

8 BB = BB +• L MU ( I , J ) *T ( J ) 

9 C ( I ) = D ( I ) * EXP ( BB ) 

WRITE! 6,530) (I, C(I), D(I), NIJ(I), 1=1, IMAX) 

WRITE! 5,540) DOSE 

400 FORMAT (21 5, El 0.4) 

401 FORMAT (7E10. 4) 

402 FORMAT ( 141 5 I 

403 FORMAT (2511) 

500 FORMAT (9H0 IC = 13/ 9H NREG = 13/ 9H IMAX = 13/ 9H DSTAR = 

1 1PE12.4) 

510 FORMAT ( 1H0/34HQREGI0N T ( I ) RHO(I) NB( I ) / ( I 7, 2F10. 3 , 1 7) ) 

530 FORMAT ( 1H0/34HC I C ( I ) D ( I ) N I J/ ( I 5 , 1 P2E i 2 . 4 , I 5 ) ) 

540 FORMAT (17HC TOTAL DOSE 1PE12.4) 

550 FORMAT ( 30H0REGI0N-J MU(I,J> ) 

570 FORMAT (16, 3X, 1P9E1 2 . 4/ ( 9X , 1P9E 1 2 . 4 ) ) 

15 RETURN 
END 
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ilBFTC INITXX 


SUBROJT INE I iMI T ( NR EG * IMAX, T, EMU, MB, C, D, A, DOS , A A, AG , G , DE S ) 
DIMENSION T ( i ) , EMU 125, 25) » NBU), CII), Dll), All), G(l), RHO(l) 
COMMON ICON, MAX, EPS, EPSD, CON, CA, WT, GG, IT, IC 

0= i.D 

DO t 0 K = 1 , MAX 
DOSA = DOS 
AAA = SQRTIAA) 

CONST = CA*Q / AAA 
IF ( D3S-DES ) 10,25,5 
5 CONST = -CONST 
10 DO 15 I =1 , MR EG 

I F ( N3 ( I ) .EQ. G) GO TO 15 
T ( I ) = T ( I ) + CONST * All) 

I F ( T ( I ) . LE , 0.0) CALL CL EAR { I , I MAX ) 

15 CONTINUE 

CALL DOZEINREG, IMAX, T, EMU, NB, C, D, A, DOS , AA, AG, G) 

I F ( ABS(DOS-DES)/DES .LE. EPSD) GO TO 25 
I F ( ( (OOS-DES)*! DOSA-OES ) ) .LE. 0.0 ) Q = 0.5 * 0 
20 CONTINUE 
K = 777 

25 CALL WEIGHT 

WRITE! 6, 30) K 

WR I TE ( 5 , 35) ! T ( I ) , I=1,NREG ) 

WRITEI6, 40) DES, DOS, WT, !D(I), 1=1, IMAX) 

RETURN 

30 FORMAT! 32H0INITIAL THICKNESSES FOUND AFTER I3.11H ITERATIONS) 

35 FORMAT! 31HGCALCULAT ED INITIAL THICKNESSES/ UP9E12.4)) 

40 FORMAT {7H0DES = 1PE12.4, 6X, 6H00S = 1PE12.4, 6X, 5HWT = 1PE12.4/ 

1 7 HOD { I ) = 1P9EI2.4/ (7X,1P9E12.4) ) 

END 


tIBFTC CLEANX 

SUBROUTINE CLEAR(J.IM) 

COMMON ICON, MAX, EPS, EPSD, CON, CA, WT, GG, IT, IC, 

1 OS T AR ( 3 ) , NREGI3), T(25,3), RH0125.3), NB!25,3), 

2 IMAX! 3 ) , D ( 2 5 , 3 ) , E MU ! 2 5, 25 , 3 ) , NIJ!25,3), D0SE13), CI25.3), 

3 A ! 2 5 , 3 ) , AA ( 3 ) , AG13), U(25,3), G!25,3) 

C THE JTH REGION HAS JUST BEEN WIPED OUT IN THE ICTH DIRECTION 

t ( j , i ; ) = o.o 

NB ( J , I C ) = 0 
DO 3 I = 1 , I M 

I F ( J .EQ. NIJ! I, IC) ) C ! 1 , 10=0.0 
5 CONTINUE 
RETURN 
END 
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c> o c? ci c~> 


SlBFTC DOSEXX 


SUBKOJTINE DOZE ( NREG , IMAX, T, EMU, NB, C, D, A, DOS , AA, AG, G) 
DIMENSION I ( i ) • EMU ( 2 5 , Z5 ) , NB(i), C(i), Dll), A(l), G ( 1 ) 

DOS = 0.0 
DO 10 1=1, IMAX 
BB - 0.0 
DO 5 J=1 , NREG 
5 BB = BB + EMU ( I, JJ* T( J) 

D C I ) = C(I) * EXP(-BB) 

10 DOS = DOS + 0 { I > 

C DOS IS TOTAL DOSE 

C AIK) = ( D-DUS E ) / ( D-X K ) 

C A A = A. A 

C AG - A . G 

AA = 0.0 

AG = 0.0 

DO 20 K = 1 , NREG 
AIK) = 0.0 

I F ( NBIK) . EQ. G) GO TO 20 
DO I 5 I = 1 , IMAX 

15 AIK) a AIK) - EMUII.K) * D(I) 

AA = AA + A I K ) ** 2 
AG = AG + A l K ) *G l K ) 

20 CONTINUE 
RETURN 
END 


$ I BFTC CYLN3D 

SUBROUTINE WEIGHT 

COMMON ICON, MAX, EPS, EPSD, CON, CA, WT, GG, IT, IC, 

1 D$TA*(3), NREG 1 3 ) , T(25,3), RHO(25,3), NB(25,3), 

2 I MAXI 3 ) , 0(25,3), EMU ( 2 5, 25 , 3 ) , NIJ(25,3>, D0SEI3), 0(25,3), 

3 A ( 2 5,3 ).j AAI3), AGO), UI25.3), GI25.3) 

DIMENSION R l 25 ) , HI25) 

GO TO (1,9, 9), ICON 

CODING FOR DOPEX SHIELDING OPTIMIZATION CODE 
SPHERICAL SHELL GEOMETRY -G.P. LAHTI -NASA-LEWIS 

FOR PLANE SLABS BOUNDED BY CONE OF HALF ANGLE THETA, 

THIS SECTIUN MAY BE USEO SIMPLY BY REPLACING 
THE TERM FORPI BY P I * I TAN ( THET A ) ) **2 
DATA FORPI / 12.56636/ 

C CALCULATE OUTER RADIUS OF EACH REGION, AND WEIGHT 

1 R 1 1 ) =T ( 1* 1 ) 

RRR= R I 1 1**3 

WT = RHO I 1,1) *RRR 
NR=NREG(1 ) 

DO : I =2, NR 

R3= RRR 

R( I )=R I 1-1 )+T( I , l) 

R RR = R I I ) ** 3 

WT= WT + RHO ( I , 1 ) • ( RRR-R3 ) 

5 CONTINUE 

WT = WT*F OR PI/3.0 
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non 


CALCULATE PARTIAL WT DERIVATIVES (DW/DXl) = GUI 

GG = j • G 

GG=:.3 

00 7 1=1, NR 

RR= 0.0 

G( I ,1) =0.0 

IFINBI I ,1 ).L 0. 0 ) GO TO 7 
DO 6 J = I , NR 
R2 = RR 

RR= R { J I **2 

6 G ( I , 1 ) =G( I»i)+RHOIJ»l)*(RR-R2) 

G { I , 1 ) =G( I , i ) » FORP I 

GG = GG + G( I ,1 ) **2 

7 CONTINUE 
RETURN 

C CYLINDER SECTION 

DATA PI / 3.14159/ 

9 NR=NRE G ( 1 ) 

R ( 1 ) =T ( 1 , 1 ) 

H ( 1 ) =T ( 1 » 2 ) 

I F ( IC0N.E0.3) GO TO 11 

DO 1C 1=2, NR 

R( I ) =R ( 1-1 ) +T ( I , 1 ) 

10 H ( I ) = H(I—1)+T(I,2) 

GO TO 14 

11 DO 12 1=2, NR 

R ( I )=R(I-1)+T( 1,1) 

12 H( I) = HI 1-1 )+T( I ,2) + T( 1,3) 

C CALCULATE WEIGHT 

14 R R=0 • 0 
WT = j • 0 

DO 15 1=1, NR 
RRR=RR 

RR= R( I )**2 * -H(I) 

WT=WT*RHO( I , 1 1 * ( RR-RRR 1 

15 CONTINUE 
WT=WT*P I 

CALCULATE PARTIAL DERIVATIVES IN THE RADIAL (FIRST) DIRECTION 
G( I ,M) = D-W/D-Tl I ,M) 

G G = ' 7 • S 
GG= 3 • 0 

DO 30 1=1, NR 
G( I , 1) =0.0 
RR=0. D 

I F ( N B ( I ,1 ) .EQ. 0) GO TO 30 
DO 25 J=I ,NR 
RRR=RR 

RR=R(J)*H( J) 

25 G(I,1)=G(I,1)+ RHO ( J , 1 ) * ( RR— RRR ) 

G ( 1,1) =G( I ,i )*2.0*PI 

GG=GG* G ( I , 1 ) **2 
30 CONTINUE 

C CALCULATE PARTIAL DERIVATIVES IN THE AXIAL DIRECTION(S) 

00 45 M=2 , ICON 
00 40 1=1, NR 
G( I , M) =0.0 
R R = 0 , D 

IFINBI I ,M) .EQ. 0) GO TO 40 
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APPENDIX B 


COMPUTER OUTPUT FROM SAMPLE PROBLEM 

TEST CYLINDER GEOMETRY 


ICON 

a 

3 

MAX 

= 

50 

EPS 

s 

1 • OOOOE-O 3 

EPSD 

= 

1. D000E-03 

CON 

= 

7.0000E-01 

CA 

= 

7.0000E-01 

IC 

3 

1 

NREG 

a 

il 

IMAX 

a 

12 

DSTAR 

a 

2.0000E+00 


REGION 

T ( I ) 

RHOU) 

NB( I ) 

I 

26.000 

9.957 

0 

2 

2.500 

8.647 

0 

3 

0.600 

16.763 

0 

4 

11.000 

9.234 

0 

5 

17.900 

0.750 

1 

6 

7.000 

19.300 

1 

7 

H.000 

0.750 

1 

8 

5.000 

19.300 

1 

9 

10.000 

0.750 

1 

10 

3.500 

19.300 

1 

11 

59.500 

0.750 

1 


REGION-J MUCI.J) 


1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 

0. 







2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 

0. 







3 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 

0. 







4 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 

C. 


0. 

0. 

0. 







5 

1.3470E-01 

1.8900E-02 

1.4700E-02 

1. 4700E-02 

2. 1200E-02 

2.4000E-02 

3.6290E-01 

2.2970E-01 

2. C180E-01 


1. 8850E-G1 

1.8370E-01 

1.7590E-01 







6 

2. 4800E-01 

7.7400E-01 

7.7200E-01 

7.7200E-01 

8.0500E-01 

7.9500E-01 

2.7190E-01 

1. 8 120E-01 

2.3860E-01 


2.5360E-01 

2. 5430E-01 

2.5770E-01 







7 

1. 3690E-01 

2 • 2800E-02 

2.0400E-02 

2.04006-02 

1. 6700E-02 

2.0600E-02 

2. 1000E-02 

2. 5000E-02 

2. 2570E-01 


1.9830E-01 

1.9100E-01 

1.7950E-01 







8 

2. 3160E-01 

7.7000E-01 

7.6600E-01 

7.6600E-01 

7.9000E-01 

7.8100E-01 

7.9800E-01 

8.27006-01 

1. 7 7S0E-01 


1.2550E-01 

2.0520E-01 

2.3680E-01 







9 

1. 2180E-01 

2.4000E-02 

2.2000E-02 

2.2000E-02 

2.0100E-02 

2.3900E-02 

1. 8700E-02 

2.45006-02 

2. 3C006-02 


2. 5000E-02 

2.2310E-01 

1 • 8080E-01 







10 

2. 3060E-01 

7.6800E-01 

7.6300E-01 

7.6300E-01 

8. 0500E-01 

7. 7600E-01 

7.8200E-01 

8. 1 200E-01 

6. ICOOE-Ol 


8. 5100E-01 

2.4070E-01 

1. 1070E-01 







11 

l. 1870E-01 

2.4400E-02 

2.2600E-02 

2.2600E-02 

1. 9600E-02 

2.5200E-02 

2. 3200E-02 

2. 8600E-02 

2. 32 COE-02 


2.9100E-02 

2. 3100E-02 

2.9700E-02 








I 

cm 

01 I > 

NIJ 

1 

2.94256*05 

2.4300E-02 

1 

2 

4.9449E*03 

3.0300E-03 

l 

3 

2. 3997E*03 

1.9600E-03 

2 

4 

2.6936E+03 

2.2000E-03 

3 

5 

3.598LE+05 

2.0400E-01 

4 

6 

1. 1522E*04 

5.0400E-03 

4 

7 

2.19736*06 

9. 2100E-02 

6 

8 

2.25216*04 

9.7400E-03 

6 

9 

9.50796*04 

9.88006-02 

8 

10 

2.0543E+04 

2.78006-02 

8 

11 

1* 1039E*05 

2.01006-01 

10 

12 

2.842SE*04 

9.47006-02 

10 


TOTAL OOSE 

7.6467E-01 



IC - 2 

NREG * XI 
I MAX * 12 

OSTAR = 2. 0000E*00 
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region 

T( II 

RHO< I ) 

NBCI» 

1 

52.000 

9.957 

0 

2 

2.500 

8.647 

0 

3 

0.600 

16.763 

0 

4 

11.000 

9.234 

0 

5 

17.900 

0.750 

1 

6 

7.000 

19.300 

1 

7 

14.000 

0.750 

1 

8 

5.000 

19.300 

1 

9 

13.000 

0.750 

1 

10 

3.500 

19.300 

1 

11 

59.500 

0.750 

1 


region-j MU ( I « J ) 


1 

0. 

0. 

0. 


0. 

0. 

0. 

2 

0. 

0. 

0. 


0. 

0. 

0. 

3 

0. 

0. 

0. 


0. 

0. 

0. 

4 

0. 

0. 

0. 


0. 

0. 

0. 

5 

1.34706-01 

1. 8900E-02 

1.4700E-02 


1. 88506-01 

1.83706-01 

1. 7590E-01 

6 

2.48006-01 

7.74006-01 

7. 7200E-01 


2. 5360E-01 

2. 54306-01 

2. 5770E-01 

7 

1. 36906-01 

2.28006-02 

2.0400E-02 


1. 9830E-01 

1.9100E-01 

1.79506-01 

8 

2. 3160E-01 

7.70006-01 

7.66006-01 


1,25506-01 

2. 05206-01 

2.36806-01 

9 

1.21806-01 

2.40006-02 

2.20006-02 


2. 50006-02 

2.2310E-01 

1.8U80E-01 

10 

2. 30606-01 

7. 6800E-01 

7.6300E-01 


8. 51006-01 

2.4070E-C1 

1. 10706-01 

n 

1.18706-01 

2.44006-C2 

2.26006-02 


2.91006-02 

2. 31006-02 

2.97006-02 


I 

cm 

0( I ) 

NIJ 

1 

2. 9<t25E+05 

2.43006-02 

1 

2 

4.9449E+03 

3.0300E-03 

1 

3 

2.3997E+03 

1.96006-03 

2 

4 

2.693SE+03 

2.2000E-03 

3 

5 

3.5981 6*05 

2.0400E-0 1 

4 

6 

1.1 5226+04 

5.0400E-03 

4 

7 

2.1 973 E+06 

9.21006-02 

6 

8 

2.25216+04 

9.74006-03 

6 

9 

9. 5079E+04 

9.8800E-02 

8 

10 

2.0543E+04 

2.78006-02 

8 

11 

1.1339E+05 

2.0100E-01 

10 

12 

2. 84256+04 

9. 4700E-02 

10 


TOTAL DOSE 

7.64676-01 



IC = 3 

nreg * 11 

IMAX » 12 

OSTAR =* 2* 30006+01 


0. 

0. 

0. 

0. 

1.4700E-02 
7.72006-01 
2. 04006-02 
7.66006-01 
2.20006-02 
7.63006-01 
2.26006-02 


0. 

0 . 

0 . 

0 . 

2.12006-02 
8.05006-01 
1 • 67006-02 
7.90006-01 
2.01006-02 
8. 0500E-01 
1.9600E-02 


0. 

0. 

0. 

0. 

2. 40006-02 
7. 9500E-01 
2.06006-02 
7.81006-01 
2.39006-02 
7.76006-01 
2.52006-02 


0. 

0. 

0. 

0. 

3.62906-01 

2.71906-01 

2.10006-02 

7.90006-01 

1.87006-02 

7.82006-01 

2.32006-02 


0. 

0. 

0. 

0. 

2.2970E-01 
1.81206-01 
2.50006-02 
8.27006-01 
2.4 500E-02 
8. 1 200E-01 
2. 86006-02 


C. 

C. 

C. 

€• 

2. 0180E-01 
2.38606-01 
2. 25706-01 
1 . 7 790E-01 
2. 30C06-02 
6. 1CC06-01 
2.22C06-02 


region 

Till 

RHOin NB(I) 









1 

52.000 

9.957 

0 









2 

2.500 

8.647 

0 









3 

3.600 

16.763 

0 









4 

11.000 

9.234 

0 









5 

17.900 

0.750 

1 









6 

7.000 

19.300 

1 









7 

14.000 

0.750 

1 









8 

5.000 

19.300 

1 









9 

13.000 

0.750 

1 









10 

3.500 

19.300 

1 









11 

59.500 

0.750 

1 









REGION-J 

MU ( I , J ) 










l 

0, 

0. 


0* 


0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 


0. 








2 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 


0. 








3 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0, 

C. 


0. 

0. 


0. 








4 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

C. 


0. 

0. 


0. 








5 

1.3470E-01 

1.89006-02 


1. 

47006-02 

1.4700E-02 

2. 1200E-02 

2.4000E-02 

3.6290E-01 

2.2970E-01 

2.01806-01 


i. 88506-01 

1.8370E-01 


1. 

75906-01 







6 

2. 4800E-01 

7. 74006-01 


7. 

72006-01 

7.72006-01 

8.0500E-01 

7.9500E-01 

2.7190E-01 

1. 81206-01 

2.38606-01 


2.53606-01 

2.5430E-01 


2. 

57706-01 







7 

1. 36906-01 

2. 28006-02 


2. 

0400E-02 

2.04006-02 

1.6700E-02 

2.06006-02 

2.1000E-02 

2. 50006-02 

2.25706-01 


1.98306-01 

1.91006-01 


1. 

7950E-01 
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8 

2* 31606-01 
1.2550E-01 

7*70006-01 

2.05206-01 

7.6600E-01 

2.3680E-01 

7* 6600E-01 

7. 9000E-01 

7.8 100E-01 

7. 9800E-01 

8.2700E-01 

1.7790E-01 

9 

1. 2180E-01 
2. 5000E-02 

2.4000E-02 

2*23106-01 

2* 2000E-02 
1*80806-01 

2.2000E-02 

2.0100E-02 

2. 3900E-02 

1. 8700E-02 

2.4500E-02 

2. 3 CC0E-02 

10 

2. 30606-01 
8.5100E-01 

7 • 6800E-01 
2.4070E-01 

7.63006-01 

1.10706-01 

7.63006-01 

8.05006-01 

7.7600E-01 

7. 82006-01 

8. 1 200E-01 

e.lOOOE-Ol 

11 

l. 16706-Oi 
2. 9100E-02 

2.4400E-G2 

2.3100E-02 

2.2600E-02 

2.97006-02 

2.2600E-02 

1. 9600E-02 

2.5200E-02 

2.320 OE-02 

2.8600E-02 

2.32C0E-C2 


I 

cm 

DI I ) 

NIJ 

1 

2.9425E+05 

2.4300E-02 

1 

2 

4. 94496+03 

3.0300E-03 

1 

3 

2.3997 E + 03 

1. 96006-03 

2 

4 

2.693SE+03 

2 . 2000E-03 

3 

5 

3.59816+05 

2 • 0400E-0 1 

4 

6 

1.1522E+04 

5.0400E-03 

4 

7 

2.19736+06 

9.2100E-02 

6 

8 

2. 25216+04 

9.7400E-03 

6 

9 

9.5L-79E+04 

9.88006-02 

8 

10 

2.05436+04 

2.7800E-02 

8 

11 

1.1039E+05 

2.01006-01 

10 

12 

2.B42SE+04 

9.4700E-02 

10 


TOTAL DOSE 

7 • 6467E-01 



THE WEIGHT OF THE STARTING CONFIGURATION IS 5.39096+07 
INITIAL THi:KN6SSES FOUND AFTER 11 ITERATIONS 
CALCULATED INITIAL THICKNESSES 

2.6000E+01 2.50006+00 6.0000E-01 l.lOOOE+Ol 1.7701E+01 6.50226+00 1.305OE+O1 4.4536E+00 9.8800E+OC 

2.8226E+O0 5.94t>8E + 01 

DES * 2 . OD 006+00 DOS = 2.0003E+00 WT * 5.01136+07 

DU)* 3.89606-02 1.15396-02 7.4045E-03 8.3112E-03 8.1703E-01 1.9633E-02 2.9957E-01 3.0622E-02 2.29376-01 

6.45596-02 3.29586-01 1.43746-01 

INITIAL THICKNESSES FOUND AFTER 11 ITERATIONS 

CALCULATED INITIAL THICKNESSES 

5.2000E+C1 2. 5000E+00 6.000C-E-01 1.10006+01 1.7701E+01 6.50226+00 1.3850E+01 4.4536E+00 9.8800E+CC 

2.8226E+J0 5. 94606+01 

OES = 2.0300E+00 DOS = 2.00036+00 WT * 4.92016+07 

DID* 3.B960E-02 1.1539E-02 7.4045E-03 8.31126-03 8.1703E-01 1.9633E-02 2.99576-01 3.0622E-02 2.29376-01 

6.45596-02 3.2958E-01 1.4374E-01 

INITIAL THICKNESSES FOUND AFTER 20 ITERATIONS 

CALCULATED INITIAL THICKNESSES 

5.20Q0E+01 2.50U0E+00 6.0000E-01 l.lOOOE+Ol 1.7467E+01 5.4B19E+00 1.3721E+01 3.3008E+00 9.7869E+CC 

!• 5166E+00 5. 9425E+01 


DES = 

2 . ODOOE+Ol 

DOS = 

2. 00026+01 

WT * 

4. 65186+07 





DC 11 = 

9 » 45 68E-02 
3. 1597E-G1 

1.7015E-01 
8. 1046E-01 

1.07586-01 
3. 0825E-01 

1.20756-01 

1.33436+01 

3.02956-01 

3.0151 E+00 

2. 932 IE-01 

1. 1201E+00 

IT = 

1 

WT 

* 4. 62366+07 

DOS - 

2.00016+00 

2. 00006+00 

2.0000E+01 



DU ) = 

4. 1448E-02 
7.33226-02 

9.8925E-03 
3. 5850E-01 

6* 3319E-03 
1.52086-01 

7. 1073E-D3 

6.9410E-01 

1. 6738E-02 

3.40046-01 

3.4888E-02 

2.65 1 5E-01 

T ( I ) = 

2.6D00E+D1 
2 . 69606+00 

2.50006+00 

5.9251E+01 

6. OOOUE-Ol 

l. 10006+01 

1. 7423E+01 

6. 8694E+00 

1 • 3573E+01 

4. 4374E+00 

9.7240 E+00 

DU ) = 

3.9S82E-02 
6. 59256-02 

1.10956-02 

3.38066-01 

7. 11446-03 
1.46396-01 

7.9855E-03 

7.8369E-01 

1.8851E-02 

3. 1001E-0 1 

3. 1674E-02 

2. 2853E-01 

tin* 

5 .2300E+01 
2.7926E+0G 

2.50006+00 
5. 9413E+01 

6.00006-01 

1. 10006+01 

1.7621E+01 

6.5937E+00 

1 . 3772E+01 

4.4498E+00 

9 • 636 C E+00 

DC I > = 

9.9193E-02 
3 • 49S9E-01 

1.64266-01 

8.68256-01 

1.0370E-01 

3.24246-01 

1.1640E-01 

1 . 2850E + 01 

2 • 9202E-01 

3.2723E+00 

3. 1576E-01 

1.2439E+00 

T( I 1 = 

5.23D0E+01 

1.43686+00 

2.5000E+00 

5.93616+01 

6. OOOOE-Ol 

1 • 1000E+01 

1. 7301E+01 

5.6021E+00 

1. 35676+01 

3. 3193 E+00 

9. 68566+00 

IT = 

2 

WT 

= 4.6016E+07 

DOS = 

2.00006+00 

2. QOOOE+OO 

2. 0000E+01 



Dill* 

4. 3952E-02 
8 .01906-02 

8.51636-03 

3.8468E-01 

5.43916-03 

1.60396-01 

6.1051E-03 

5.9166E-01 

1.43346-02 

3.7422E-01 

3.8691E-02 

2.S187E-01 
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T 1 1 > = 

2.5300E+Oi 
Z, 6299E+CU 

2.5000E+00 

5.8987E+01 

6. 0000E-01 

1.10006+01 

1. 7201E+01 

7# 2050E+00 

1.3310E+01 

4. 3866E+00 

9.*7e5E+00 

Dtl ) = 

4.0542E-02 

6.9534E-02 

1.0586E-02 

3.4799E-01 

6.7828E-03 
X. 4958E-01 

7.6134E-03 

7.4561E-01 

1. 7958E-02 

3.2212E-01 

3.2879E-02 

2.488CE-01 

tid* 

5.2000E+01 
2 . 76 31 E ♦00 

2.5000E+00 

5.93456+01 

6. OOOOE-Ol 

1.10006+01 

1.75276+01 

6. 6935E+00 

1.3681E+01 

4.4435E+00 

9.78646+00 

om = 

1 .05046-01 
3 • 9332E-01 

1.56846-01 
9. 4320E-OX 

9. 8847E-02 
3.4403E-OX 

1.10956-01 

1.22326+01 

2.7832E-01 

3* 5889E+00 

3.4277E-01 

1.4050E+00 

TIM* 

5 • 2300E+01 
1. 34796+00 

2.5000E+00 

5.9280E+01 

6.0000E-01 

1.1000E+01 

1.7100E+01 

5. 7416E+00 

1.3379EV01 

3.34466+00 

9.5623E+00 

IT * 

3 

WT 

= 4.5834E+07 

DOS = 

2.00006+00 

2. OOOOE+OO 

2.0000E+01 



Oil)® 

4.5572E-02 
8.23 95E-02 

7. 56866-03 
4.0645E-01 

4.8252E-03 

1.68806-01 

5# 4161 E-03 

5.2074E-01 

1.2686E-02 

3.9952E-0X 

4. 1826E-02 

3.035 IE— 01 

T( I ) = 

2 » 6000E+01 
2 • 63 58E+00 

2.5000E+00 

5.8670E+01 

6. 00006-01 

1.1000E+01 

1.70666+01 

7. 4748E+00 

1.3076E»01 

4. 28696+00 

9. 4549E+00 

om* 

4. 1547E-02 
7.22I1E-02 

1.00246-02 
3. 5936E-01 

6.41666-03 

1.5336E-01 

7. 20236-03 

7.03566-01 

1.6975E-02 

3.3551E-01 

3.41896-02 

2.59636-01 

Tm = 

5 .20006+01 
2 • 73 77E+00 

2.5000E+00 

5.9262E+01 

6. OOOOE-Ol 

1.10006+01 

1. 7419E+01 

6. 8144E+00 

1.3S80E+01 

4.43336+00 

9. 73146+00 

D ( I ) = 

I.1235E-01 
4* 4B36E-01 

1.4775E-01 

1.0398E+00 

9.2930E-02 
3. 7122E-01 

1.04316-01 

1.14786+01 

2. 6163E-01 

3# 9609E+00 

3.73306-01 

1 . 6098E+00 

Tin* 

5. 20006+01 
I.2549E+0O 

2.5000E+00 

5.9182E+01 

6. OOOOE-OX 

1.1000E+01 

1.6864E+01 

5.8969E+00 

X.3I55E+0I 

3. 3789E+00 

9.41 37E+00 

IT * 

4 

WT 

* 4.5675E+07 

00S = 

2. OOOOE+OO 

2. OOOOE+OO 

2.0000E+01 



D(I )* 

4# 93306-02 
8 .02066—02 

7.0540E-03 
4.2281 E— 01 

4. 4911E-03 
1.76956-01 

5. 041 IE-03 

4.0121E-O1 

1. 17896-02 

9.1551E-01 

4.4243E-02 

3.C136E-01 

T 1 1 ) = 

2.60006+01 
2 .69766+00 

2. 5000E+00 
5.8309E+QX 

6. OOOOE-Ol 

1. 1000E+01 

1.70256+01 

7.66206+00 

X.2878E.0X 

4. 1497E+00 

9. 3641E+00 

dim* 

4.2S50E-02 

7.4586E-02 

9.4576 E-03 
3.7148E-01 

6. 0487E-03 
1.5760E-01 

6. 78936-03 

6.6132E-01 

1.59886-02 

3.4893E-01 

3.5482E-02 

2.6968E-01 

Tm* 

5 » 2300E + 01 
2 . 72 I5E+00 

2.5000E+00 

5.9167E+01 

6. OOOOE-Ol 

1.1000E+01 

1. 73076+01 

6.9314E+00 

I.3474E+0X 

4.4181E+00 

9.6736E+00 

did* 

X.2394E-0X 

5.X193E-0X 

1.37626-01 
1 • X 574E+00 

8.6371E-02 

4.0330E-01 

9.6947E-02 

1.06416+01 

2.4314E-01 

8.3496E+00 

4.0346E-01 

1 • 6489E+00 

Tin* 

5.2300E+OX 
X. IS99E+0U 

2.5000E+00 
5. 9069E+01 

6. OOOOE-Ol 

1.10006+01 

1.66096+01 

6. 0530E+00 

1 • 2905E+01 

3.4216E+00 

9. 2438E+00 

IT * 

5 

WT 

= 4. 5529E+07 

00S * 

2. OOOOE+OO 

2. OOOOE+OO 

2.0000E+OX 



D(I ) = 

5 . 22 01E-02 
7.63626-C2 

6.8152E-03 

4.3452E-01 

4.3346E-03 
1 . 8428E-01 

4.86546-03 

4.6154E-01 

1.1368E-02 

4.2903E-01 

4.6079E-02 

2.9310E-01 

TCI ) = 

2.6000E+01 
2 • 7301E+00 

2.5000E+00 

5.7916E+01 

6. OOOOE-Ol 

1. 10006+01 

1.7055E+01 

7. 7880E+00 

1.2706E+01 

4.0041E+00 

9. 2059E+00 

did* 

4.3300E-02 
7 . 6375E-02 

8.9340E-03 

3.8363E-01 

5.7087E-03 
1 • 6209E-01 

6.4078E-03 

6.2226E-01 

1.5078E-02 

3*6X24E-0X 

3.66526-02 

2.7783E-01 

Tin* 

5 • 23 OOE+01 
2* 7170E+00 

2. 5000E+00 
5.9063E+01 

6. OOOOE-OX 

1* 1000E+01 

1.7200E+01 

7. 0408E+00 

X.337XE+0X 

4. 3972E+00 

9. 6 15 CE+00 

did* 

1.3055E-01 

5.7378E-01 

1.2735E-01 

1.2933E+00 

7.9758E-02 
4. 4059E-01 

8.9525E-02 

9.7958E+00 

2. 2452E-01 

4.7074E*00 

4. 2898E-01 

2 • 1036E+00 

Tin* 

5.2300E+CX 
1 . 1036E+00 

2.5000E+00 

5.8945E+0X 

6. COOOE-Ol 

1. 1000E+01 

1. 6354E+01 

6. 1944E+00 

X.2643E+0X 

3. 4699E+00 

9.C587E+00 

IT * 

29 

WT 

* 4.3 547E+07 

DOS * 

2. OOOOE+OO 

2. OOOOE+OO 

2.OOOOE+0X 



did* 

1.99326-01 

6.12B9E-G2 

7.2807E-03 

3.7302E-01 

4.5505E-03 

2.0604E-01 

5. 10776-03 

4.4022E-01 

1.18466-02 

3*8277 E-OX 

7.1793E-02 

2 . 5675E-01 

Tin* 

2. 6300E+0 L 

2.5000E+00 

6. QOOOE-Oi 

l. 1000E+0L 

1.98086+01 

9* 1816E+00 

9*8857E*00 

2.35706+00 

1.C864E+01 


3 . 33 85E+00 4.4377E+01 
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D( I ) = 

6.8322E-02 
7 ; 23 22E-02 

7.06036-03 

*.*717E-01 

4.47*66-03 

1.97266-01 

5 ■ 0225E-03 

4.7185E-01 

1. 1823E-02 

3.99306-01 

4. 3264E-02 

2.7163 E- 01 

T ( I 1 = 

5 » 2000E+01 
2.9480E*G0 

2.50006*00 
5. 5394E+01 

6.0000E-01 

1.10006*01 

1.70846*01 

7.62966*00 

1.28266*01 

4.02606*00 

9 • 29 1 5E*00 

0(11 = 

4.3410E-01 

6.82336-01 

6.8282E-02 

*.887*6*00 

4. 1856E-02 
1.6638E*00 

4.69826-02 

4.82216*00 

1.20036-01 

4.1679E*00 

3.02096-01 

2.76316*00 

r (i ) = 

5.2D0OE+O1 

2.2D62E+00 

2. 5000E*00 
5. *815E*01 

6. 0000E-01 

1. 1000E+01 

1. 4612E*01 

6.08996*00 

9.7294E*00 

3. 6326E + 00 

6 • 23 58E*00 

IT = 

30 

HT 

= 4.3*966*07 

DOS = 

2 • OQO0E*O0 

2. 0000E*00 

2.0000E*01 



DU 1 = 

2.0945E-01 

6.1571E-02 

7.3122 E-03 
3.6*736-01 

4.56726-03 
2. 02916-01 

5.12656-03 

4.40506-01 

1. 1887E-02 

3.61216-01 

7. 3260E-02 

2 • 5748E-0 1 

Til ) = 

2.63006*01 

3.3240E*00 

2.50006*00 
*• 3763 E*0i 

6. 000QE-01 

1.10006*01 

1.99206*01 

9.23946*00 

9.84516*00 

2. 32286*00 

1 • C952E*01 

DU ) = 

6.9740E-02 
7. 33946-02 

7.10166-03 

*.*3266-01 

4.50006-03 

1.95216-01 

5.0510E-03 

4.74526-01 

1.18936-02 

3.9897E-01 

4. 3419E-02 

2.73246-01 

T(I >* 

5.23306*01 
2 • 92 23E*00 

2.50006*00 

5.52016*01 

6. 00006-01 

1.10006*01 

1.71076*01 

7.6327E*00 

1.28726*01 

4. 0449E*00 

9 • 299CE+00 

DU ) = 

4.4145E-01 
6. 85886-01 

6.86136-02 

4.85856*00 

4.20506-02 

1.653*6*00 

4.71996-02 

4.8429E*Q0 

i. 20626-01 

4. 1640 E*00 

3 • 0307E-01 

2 • 77 14E*00 

T 1 1 1 = 

5 • 2300E+01 
2.19006^-00 

2.50006*00 

5.46086*01 

6.00006-01 

1.10006*01 

1.46336*01 

6.0932E+00 

9. 7640E+00 

3.64416*00 

6 . 240 4E*00 

IT = 

31 

HT 

= 4.3**86*07 

DOS = 

2. 0000E*00 

2.00006*00 

2. OOOOE+O 1 



D ( I ) = 

2 . 1983E-01 
6.19736-02 

7. 3434E-03 
3.5623E-01 

4. 5838E-03 
1.99536-01 

5. 14516-03 

4.4079E-01 

1. 19276-02 

3. 5977E-01 

7.4793E-02 

2.58196-01 

Til I* 

2.6300E*01 

3.3367E*00 

2.5000E*00 
*. 3151 E*01 

6.00006-01 

1 . 1000E*01 

2. 0032E*01 

9. 29836*00 

9.8120E+00 

2.2903E+00 

1 . 1038E+01 

DID* 

7 * 0563E-02 
7.39186-02 

7.1437E-03 

4.39116-01 

4. 5259E-03 
1 • 9301E-01 

5.08016-03 

4. 77276-01 

1.1954E-02 

3.98816-01 

4.35886-02 

2. 7492E-01 

TUI* 

5. 20006*01 
2.89*56*00 

2* 5000E+00 
5.50056*01 

6.00006-01 

1.10006*01 

1.71306*01 

7.63596*00 

1.2925E*0L 

4.06576*00 

9 • 3056 E+CO 

DU I* 

*.*3706-01 

6.9254E-01 

6.89706-02 

4.82336*00 

4.2260E-02 

1.6401E*00 

4.74346-02 

4.86596*00 

1.21266-01 

4.16276*00 

3. 04296-01 

2.78266*00 

Till* 

5.2000E*01 
2. 15986*00 

2.50006*00 

5.43966*01 

6.00006-01 

1.10006*01 

1.46546*01 

6. 097 2E*00 

9.80556*00 

3.65836*00 

6.2461 E*C 0 

IT * 

32 

HT 

*> *.3*01E*07 

DOS = 

2.0000E*00 

2.00006*00 

2.00006*01 



DID* 

2 . 3341E-01 
6. 21 93E-02 

7.37*26-03 

3.47566-01 

4. 60026-03 
1.9591E-01 

5. 16356-03 

4. 4110E-01 

1. 19676-02 

3.58456-01 

7. 63976-02 

2. 5888E-01 

TUI* 

2.63306*01 

3.20656*00 

2.50006*00 
*• 25*5E*01 

6.00006-01 

1. 10006*01 

2.01456*01 

9.35846*00 

9.78636*00 

2.25916*00 

1. il23E*01 

DU I* 

7 • 1592E-02 
7. *7806-02 

7.1865E-03 
*• 3473E-01 

4.55236-03 

1.9066E-01 

5. 10976-03 

4.80106-01 

1.20376-02 

3.98826-01 

4.37716-02 

2. 7665E-01 

TUI* 

5. 20306*01 
2.86*76*00 

2.50006*00 

5.48046*01 

6.00006-01 

1.10006*01 

1.71516*01 

7.63926*00 

1.29846*01 

4.08836*00 

9.311CE*00 

DU 1 = 

*.55916-01 

6.991*6-01 

6.93*96-02 

4.78296*00 

4.24836-02 

1.62446*00 

4. 76856-02 

4.89086*00 

1.21936-01 

4.1636 6*00 

3. 057 CE- 01 

2.79616400 

TUI* 

5.23006*01 

2.1*616*00 

2.50006*00 

5.41806*01 

6.00006-01 

1.10006*01 

1.46756*01 

6.10186*00 

9.8540E*00 

3.67496*00 

6 . 2 52 2E*00 

IT * 

33 

HT 

* 4.33576*07 

DOS - 

2.00006*00 

2.00006*00 

2.00006*01 



DU 1“ 

2. *1186-01 
6. 2528E-02 

7.4048E-03 

3.38746-01 

4.61656-03 

1.92066-01 

5. 1818E-03 

4.41436-01 

1.20056-02 

3.57256-01 

7.8078E-02 

2 • 5953E-01 

TUI* 

2 . 63 30E+01 
3.26366*00 

2.5000E*00 
4. 19436*01 

6.00006-01 

1.10006*01 

2. 0259E*01 

9.4198E*00 

9.76796*00 

2.22926*00 

1 .12046*01 

D(I I* 

7.25326-02 

7.57006-02 

7.23006-03 
4. 3014E-01 

4.57906-03 

1.88186-01 

5.13976-03 

4.82996-01 

1.21126-02 

3.99026-01 

4.39716-02 

2 • 704 IE-01 


29 



T(I) = 

5.23DQE+OI 
2. 8330E+00 

2 • 5000E+00 
5.4599E+01 

6.0000E-01 

1 • 1000E+01 

0(11 = 

4.6315E-0I 
7 » 05 54E-01 

6.9745E-02 
4. 7382 E^OO 

4.2716E-02 
1. 6066E+00 

4.7947E-02 

1111 = 

5 • 2D 30E + 01 
2.1L92E+00 

2.5000E+00 

5.3957E+01 

6.0000E-01 

1. lOOOE+Ol 

IT = 

34 WT 

= 4.3314E+07 

DOS = 

D(I ( = 

2. 5237E-01 
6. 2876E-02 

7.4352E-03 

3.29806-01 

4.6328E-03 
1. 8800E-01 

5.2001E-03 

Ttl> = 

2.to;>oE*oi 
3 . 2378E+00 

2.50006400 
4. 13476+01 

fe.OOOOE-Ol 

l.lOOOEvOi 

0(11 = 

7.3V84E-02 

7.6S50E-02 

7 • 2743E-03 
4.2534E-0I 

4.6062E-03 
1* 8557E-01 

5.1703E-03 

nn= 

5 » 23 OOE+Ol 
2 . 7993E+00 

2. 50006*00 
5.43886+01 

6.0000E-01 

1 • 1000E+Q1 

0(11 = 

4 • 73 %5E-0i 
7.1451E-01 

7.0157E-02 

4.68976+00 

4.2959E-02 
1* 5871E+00 

4.8220E-02 

T ( I I s 

5 . 2D DOE+Oi 
2.0895E+00 

2. 5000E+00 
5.37296+01 

6.0000E-01 

1 • 1000E+01 


1.7171E+01 

7.6426E+00 

1 • 3050E+01 

4. 11266*00 

9 • 21 49E+00 

4.9173E+00 

1 .22646-01 

4.1665E+00 

3.0729E-01 

2.81 12E+00 

1.4696E+01 

6. 1069E+00 

9. 90956*00 

3 • 6939E+00 

6 • 258 2 E+00 

! • OOOOE+OO 

2. OOOOE+OO 

2.0000E+01 



4.4178E-01 

1 • 2043E-02 

3.5616E-01 

7 • 9845E-02 

2.6014E-01 

2.0374EV01 

9. 4826E+00 

9.7567E*00 

2 • 2004E*00 

1 • 12B4E*01 

4.8595E-01 

1.21886-02 

3.99426-01 

4.4187E-02 

2.aoieE-oi 

1.7190E4-01 

7.6461E+00 

1. 31236*01 

4. 13876*00 

9. 3172E*00 

4.9452E+00 

l* 2338E-01 

4.17156*00 

3.0904E-01 

2 • £2 7 5E*0G 

1.4717E+01 

6* 1124E + 00 

9.9723E*00 

3. 71506*00 

6 . 2634E*00 
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